Pathological claustrum activity drives aberrant cognitive network processing in human chronic pain.

Aberrant cognitive network activity and cognitive deficits are established features of chronic pain. However, the nature of cognitive network alterations associated with chronic pain and their underlying mechanisms require elucidation. Here, we report that the claustrum, a subcortical nucleus implicated in cognitive network modulation, is activated by acute painful stimulation and pain-predictive cues in healthy participants. Moreover, we discover pathological activity of the claustrum and a region near the posterior inferior frontal sulcus of the right dorsolateral prefrontal cortex (piDLPFC) in migraine patients during acute pain and cognitive task performance. Dynamic causal modeling suggests a directional influence of the claustrum on activity in this piDLPFC region, and diffusion weighted imaging verifies their structural connectivity. These findings advance understanding of claustrum function during acute pain and provide evidence of a possible circuit mechanism driving cognitive impairments in chronic pain.

Pathological claustrum activity drives aberrant cognitive network processing in human chronic pain.

Aberrant cognitive network activity and cognitive deficits are established features of chronic pain. However, the nature of cognitive network alterations associated with chronic pain and their underlying mechanisms require elucidation. Here, we report that the claustrum, a subcortical nucleus implicated in cognitive network modulation, is activated by acute painful stimulation and pain-predictive cues in healthy participants. Moreover, we discover pathological activity of the claustrum and a lateral aspect of the right dorsolateral prefrontal cortex (latDLPFC) in migraine patients. Dynamic causal modeling suggests a directional influence of the claustrum on activity in this latDLPFC region, and diffusion weighted imaging (DWI) verifies their structural connectivity. These findings advance understanding of claustrum function during acute pain and provide evidence of a possible circuit mechanism driving cognitive impairments in chronic pain.

Regression models for partially localized fMRI connectivity analyses.

Background: Brain functional connectivity analysis of resting-state functional magnetic resonance imaging (fMRI) data is typically performed in a standardized template space assuming consistency of connections across subjects. Analysis methods can come in the form of one-edge-at-a-time analyses or dimension reduction/decomposition methods. Common to these approaches is an assumption that brain regions are functionally aligned across subjects; however, it is known that this functional alignment assumption is often violated. Methods: In this paper, we use subject-level regression models to explain intra-subject variability in connectivity. Covariates can include factors such as geographic distance between two pairs of brain regions, whether the two regions are symmetrically opposite (homotopic), and whether the two regions are members of the same functional network. Additionally, a covariate for each brain region can be included, to account for the possibility that some regions have consistently higher or lower connectivity. This style of analysis allows us to characterize the fraction of variation explained by each type of covariate. Additionally, comparisons across subjects can then be made using the fitted connectivity regression models, offering a more parsimonious alternative to edge-at-a-time approaches. Results: We apply our approach to Human Connectome Project data on 268 regions of interest (ROIs), grouped into eight functional networks. We find that a high proportion of variation is explained by region covariates and network membership covariates, while geographic distance and homotopy have high relative importance after adjusting for the number of predictors. We also find that the degree of data repeatability using our connectivity regression model-which uses only partial location information about pairs of ROI's-is comparably as high as the repeatability obtained using full location information. Discussion: While our analysis uses data that have been transformed into a common template-space, we also envision the method being useful in multi-atlas registration settings, where subject data remains in its own geometry and templates are warped instead. These results suggest the tantalizing possibility that fMRI connectivity analysis can be performed in subject-space, using less aggressive registration, such as simple affine transformations, multi-atlas subject-space registration, or perhaps even no registration whatsoever.

Sliding windows analysis can undo the effects of preprocessing when applied to fMRI data.

Resting-state fMRI (rs-fMRI) data is used to study the intrinsic functional connectivity (FC) in the human brain. Recently, interest has focused on studying the temporal dynamics of FC on short timescales, ranging from seconds to minutes. These studies of time-varying FC (TVFC) have enabled the classification of whole-brain dynamic FC profiles into distinct "brain states", defined as recurring whole-brain connectivity profiles reliably observed across subjects and sessions. The analysis of rs-fMRI data is complicated by the fact that the measured BOLD signal consists of changes induced by neuronal activation, as well as non-neuronal nuisance fluctuations that should be removed prior to further analysis. Thus, the data undergoes significant preprocesing prior to analysis. In previous work [19], we illustrated the potential pitfalls involved with using modular preprocessing pipelines, showing how later preprocessing steps can reintroduce signal variation previously removed from the data. Here we show that the problem runs deeper, and that certain statistical analysis techniques can similarly interact with preprocessing and reintroduce previously removed signal variation. One such technique is the popular sliding window analysis, used to compute TVFC. In this paper, we discuss the problem both theoretically and empirically in application to test-retest rs-fMRI data. Importantly, we show that we are able to obtain essentially the same brain states and state transitions when analyzing motion induced signal as we do when analyzing the preprocessed but windowed data. Our results cast doubt on whether the estimated brain states obtained using sliding window analysis are neuronal in nature, or simply reflect non-neuronal nuisance signal variation (e.g., motion). We highlight the need to critically revisit previous work on rs-fMRI data work that may not have adequately controlled for these types of effects.

Placebo treatment affects brain systems related to affective and cognitive processes, but not nociceptive pain.

Placebo analgesia is a replicable and well-studied phenomenon, yet it remains unclear to what degree it includes modulation of nociceptive processes. Some studies find effects consistent with nociceptive effects, but meta-analyses show that these effects are often small. We analyzed placebo analgesia in a large fMRI study (N = 392), including placebo effects on brain responses to noxious stimuli. Placebo treatment caused robust analgesia in both conditioned thermal and unconditioned mechanical pain. Placebo did not decrease fMRI activity in nociceptive pain regions, including the Neurologic Pain Signature (NPS) and pre-registered spinothalamic pathway regions, with strong support from Bayes Factor analyses. However, placebo treatment affected activity in pre-registered analyses of a second neuromarker, the Stimulus Intensity Independent Pain Signature (SIIPS), and several associated a priori brain regions related to motivation and value, in both thermal and mechanical pain. Individual differences in behavioral analgesia were correlated with neural changes in both thermal and mechanical pain. Our results indicate that processes related to affective and cognitive aspects of pain primarily drive placebo analgesia.

Improved fMRI-based pain prediction using Bayesian group-wise functional registration.

In recent years, the field of neuroimaging has undergone a paradigm shift, moving away from the traditional brain mapping approach towards the development of integrated, multivariate brain models that can predict categories of mental events. However, large interindividual differences in both brain anatomy and functional localization after standard anatomical alignment remain a major limitation in performing this type of analysis, as it leads to feature misalignment across subjects in subsequent predictive models. This article addresses this problem by developing and validating a new computational technique for reducing misalignment across individuals in functional brain systems by spatially transforming each subject's functional data to a common latent template map. Our proposed Bayesian functional group-wise registration approach allows us to assess differences in brain function across subjects and individual differences in activation topology. We achieve the probabilistic registration with inverse-consistency by utilizing the generalized Bayes framework with a loss function for the symmetric group-wise registration. It models the latent template with a Gaussian process, which helps capture spatial features in the template, producing a more precise estimation. We evaluate the method in simulation studies and apply it to data from an fMRI study of thermal pain, with the goal of using functional brain activity to predict physical pain. We find that the proposed approach allows for improved prediction of reported pain scores over conventional approaches. Received on 2 January 2017. Editorial decision on 8 June 2021.

Evaluating the prediction performance of objective physical activity measures for incident Parkinson's disease in the UK Biobank.

BACKGROUND: Parkinson's disease (PD) is the fastest-growing neurological condition with over 10 million cases worldwide. While age and sex are known predictors of incident PD, there is a need to identify other predictors. This study compares the prediction performance of accelerometry-derived physical activity (PA) measures and traditional risk factors for incident PD in the UK Biobank. METHODS: The study population consisted of 92,352 UK Biobank participants without PD at baseline (43.8% male, median age 63 years with interquartile range 43-69). 245 participants were diagnosed with PD by April 1, 2021 (586,604 person-years of follow-up). The incident PD prediction performances of 10 traditional predictors and 8 objective PA measures were compared using single- and multi-variable Cox models. Prediction performance was assessed using a novel, stable statistic: the repeated cross-validated concordance (rcvC). Sensitivity analyses were conducted where PD cases diagnosed within the first six months, one year, and two years were deleted. RESULTS: Single-predictor Cox regression models indicated that all PA measures were statistically significant (p-values < 0.0001). The highest-performing individual predictors were total acceleration (TA) (rcvC = 0.813) among PA measures, and age (rcvC = 0.757) among traditional predictors. The two-step forward-selection process produced a model containing age, sex, and TA (rcvC = 0.851). Adding TA to the model increased the rcvC by 9.8% (p-value < 0.0001). Results were largely unchanged in sensitivity analyses. CONCLUSIONS: Objective PA summaries have better single-predictor model performance than known risk factors and increase the prediction performance substantially when added to models with age and sex.

Brain mediators of negative affect-induced physical symptom reporting in patients with functional somatic syndromes.

Functional somatic syndromes (FSS) include fibromyalgia, irritable bowel syndrome (IBS), and others. In FSS patients, merely viewing negative affective pictures can elicit increased physical symptoms. Our aim was to investigate the neural mechanisms underlying such negative affect-induced physical symptoms in FSS patients. Thirty patients with fibromyalgia and/or IBS and 30 healthy controls (all women) watched neutral, positive and negative affective picture blocks during functional MRI scanning and rated negative affect and physical symptoms after every block. We compared brain-wide activation during negative versus neutral picture viewing in FSS patients versus controls using robust general linear model analysis. Further, we compared neurologic pain signature (NPS), stimulus intensity-independent pain signature (SIIPS) and picture-induced negative emotion signature (PINES) responses to the negative versus neutral affect contrast and investigated whether they mediated between-group differences in affective picture-induced physical symptom reporting. More physical symptoms were reported after viewing negative compared to neutral pictures, and this effect was larger in patients than controls (p = 0.025). Accordingly, patients showed stronger activation in somatosensory regions during negative versus neutral picture viewing. NPS, but not SIIPS nor PINES, responses were higher in patients than controls during negative versus neutral pictures (p = 0.026). These differential NPS responses partially mediated between-group differences in physical symptoms. In conclusion, picture-induced negative affect elicits physical symptoms in FSS patients as a result of activation of somatosensory and nociceptive brain patterns, supporting the idea that affect-driven alterations in processing of somatic signals is a critical mechanism underlying FSS.

Quantifying the Association between Objectively Measured Physical Activity and Multiple Sclerosis in the UK Biobank.

INTRODUCTION: Objectively measured physical activity (PA) data were collected in the accelerometry substudy of the UK Biobank. UK Biobank also contains information about multiple sclerosis (MS) diagnosis at the time of and after PA collection. This study aimed to 1) quantify the difference in PA between prevalent MS cases and matched healthy controls, and 2) evaluate the predictive performance of objective PA measures for incident MS cases. METHODS: The first analysis compared eight accelerometer-derived PA summaries between MS patients ( N = 316) and matched controls (30 controls for each MS case). The second analysis focused on predicting time to MS diagnosis among participants who were not diagnosed with MS. A total of 19 predictors including eight measures of objective PA were compared using Cox proportional hazards models (number of events = 47; 585,900 person-years of follow-up). RESULTS: In the prevalent MS study, the difference between MS cases and matched controls was statistically significant for all PA summaries ( P < 0.001). In the incident MS study, the most predictive variable of progression to MS in univariate Cox regression models was lower age ( C = 0.604), and the most predictive PA variable was lower relative amplitude (RA, C = 0.594). A two-stage forward selection using Cox regression resulted in a model with concordance C = 0.693 and four predictors: age ( P = 0.015), stroke ( P = 0.009), Townsend deprivation index ( P = 0.874), and RA ( P = 0.004). A model including age, stroke, and RA had a concordance of C = 0.691. CONCLUSIONS: Objective PA summaries were significantly different and consistent with lower activity among study participants who had MS at the time of the accelerometry study. Among individuals who did not have MS, younger age, stroke history, and lower RA were significantly associated with a higher risk of a future MS diagnosis.

A dorsomedial prefrontal cortex-based dynamic functional connectivity model of rumination.

Rumination is a cognitive style characterized by repetitive thoughts about one's negative internal states and is a common symptom of depression. Previous studies have linked trait rumination to alterations in the default mode network, but predictive brain markers of rumination are lacking. Here, we adopt a predictive modeling approach to develop a neuroimaging marker of rumination based on the variance of dynamic resting-state functional connectivity and test it across 5 diverse subclinical and clinical samples (total n = 288). A whole-brain marker based on dynamic connectivity with the dorsomedial prefrontal cortex (dmPFC) emerges as generalizable across the subclinical datasets. A refined marker consisting of the most important features from a virtual lesion analysis further predicts depression scores of adults with major depressive disorder (n = 35). This study highlights the role of the dmPFC in trait rumination and provides a dynamic functional connectivity marker for rumination.

Predicting chronic postsurgical pain: current evidence and a novel program to develop predictive biomarker signatures.

ABSTRACT: Chronic pain affects more than 50 million Americans. Treatments remain inadequate, in large part, because the pathophysiological mechanisms underlying the development of chronic pain remain poorly understood. Pain biomarkers could potentially identify and measure biological pathways and phenotypical expressions that are altered by pain, provide insight into biological treatment targets, and help identify at-risk patients who might benefit from early intervention. Biomarkers are used to diagnose, track, and treat other diseases, but no validated clinical biomarkers exist yet for chronic pain. To address this problem, the National Institutes of Health Common Fund launched the Acute to Chronic Pain Signatures (A2CPS) program to evaluate candidate biomarkers, develop them into biosignatures, and discover novel biomarkers for chronification of pain after surgery. This article discusses candidate biomarkers identified by A2CPS for evaluation, including genomic, proteomic, metabolomic, lipidomic, neuroimaging, psychophysical, psychological, and behavioral measures. Acute to Chronic Pain Signatures will provide the most comprehensive investigation of biomarkers for the transition to chronic postsurgical pain undertaken to date. Data and analytic resources generatedby A2CPS will be shared with the scientific community in hopes that other investigators will extract valuable insights beyond A2CPS's initial findings. This article will review the identified biomarkers and rationale for including them, the current state of the science on biomarkers of the transition from acute to chronic pain, gaps in the literature, and how A2CPS will address these gaps.

Regression models for partially localized fMRI connectivity analyses.

Brain functional connectivity analysis of resting-state functional magnetic resonance imaging (fMRI) data is typically performed in a standardized template space assuming consistency of connections across subjects. This can come in the form of one-edge-at-a-time analyses or dimension reduction/decomposition methods. Common to these approaches is the assumption of complete localization (or spatial alignment) of brain regions across subjects. Alternative approaches completely eschew localization assumptions by treating connections as statistically exchangeable (for example, using the density of connectivity between nodes). Yet other approaches, such as hyperalignment, attempt to align subjects on function as well as structure, thereby achieving a different sort of template-based localization. In this paper, we propose the use of simple regression models to characterize connectivity. To that end, we build regression models on subject-level Fisher transformed regional connection matrices using geographic distance, homotopic distance, network labels, and region indicators as covariates to explain variation in connections. While we perform our analysis in template-space in this paper, we envision the method being useful in multi-atlas registration settings, where subject data remains in its own geometry and templates are warped instead. A byproduct of this style of analysis is the ability to characterize the fraction of variation in subject-level connections explained by each type of covariate. Using Human Connectome Project data, we found that network labels and regional characteristics contribute far more than geographic or homotopic relationships (considered non-parametrically). In addition, visual regions had the highest explanatory power (i.e., largest regression coefficients). We also considered subject repeatability and found that the degree of repeatability seen in fully localized models is largely recovered using our proposed subject-level regression models. Further, even fully exchangeable models retain a sizeable amount of repeatability information, despite discarding all localization information. These results suggest the tantalizing possibility that fMRI connectivity analysis can be performed in subject-space, using less aggressive registration, such as simple affine transformations, multi-atlas subject-space registration, or perhaps even no registration whatsoever.

Genetic risk shared across 24 chronic pain conditions: identification and characterization with genomic structural equation modeling.

ABSTRACT: Chronic pain conditions frequently co-occur, suggesting common risks and paths to prevention and treatment. Previous studies have reported genetic correlations among specific groups of pain conditions and reported genetic risk for within-individual multisite pain counts (≤7). Here, we identified genetic risk for multiple distinct pain disorders across individuals using 24 chronic pain conditions and genomic structural equation modeling (Genomic SEM). First, we ran individual genome-wide association studies (GWASs) on all 24 conditions in the UK Biobank ( N ≤ 436,000) and estimated their pairwise genetic correlations. Then we used these correlations to model their genetic factor structure in Genomic SEM, using both hypothesis- and data-driven exploratory approaches. A complementary network analysis enabled us to visualize these genetic relationships in an unstructured manner. Genomic SEM analysis revealed a general factor explaining most of the shared genetic variance across all pain conditions and a second, more specific factor explaining genetic covariance across musculoskeletal pain conditions. Network analysis revealed a large cluster of conditions and identified arthropathic, back, and neck pain as potential hubs for cross-condition chronic pain. Additionally, we ran GWASs on both factors extracted in Genomic SEM and annotated them functionally. Annotation identified pathways associated with organogenesis, metabolism, transcription, and DNA repair, with overrepresentation of strongly associated genes exclusively in brain tissues. Cross-reference with previous GWASs showed genetic overlap with cognition, mood, and brain structure. These results identify common genetic risks and suggest neurobiological and psychosocial mechanisms that should be targeted to prevent and treat cross-condition chronic pain.

Associations between neighborhood socioeconomic status, parental education, and executive system activation in youth.

Socioeconomic status (SES) can impact cognitive performance, including working memory (WM). As executive systems that support WM undergo functional neurodevelopment during adolescence, environmental stressors at both individual and community levels may influence cognitive outcomes. Here, we sought to examine how SES at the neighborhood and family level impacts task-related activation of the executive system during adolescence and determine whether this effect mediates the relationship between SES and WM performance. To address these questions, we studied 1,150 youths (age 8-23) that completed a fractal n-back WM task during functional magnetic resonance imaging at 3T as part of the Philadelphia Neurodevelopmental Cohort. We found that both higher neighborhood SES and parental education were associated with greater activation of the executive system to WM load, including the bilateral dorsolateral prefrontal cortex, posterior parietal cortex, and precuneus. The association of neighborhood SES remained significant when controlling for task performance, or related factors like exposure to traumatic events. Furthermore, high-dimensional multivariate mediation analysis identified distinct patterns of brain activity within the executive system that significantly mediated the relationship between measures of SES and task performance. These findings underscore the importance of multilevel environmental factors in shaping executive system function and WM in youth.

Enhanced hyperalignment via spatial prior information.

Functional alignment between subjects is an important assumption of functional magnetic resonance imaging (fMRI) group-level analysis. However, it is often violated in practice, even after alignment to a standard anatomical template. Hyperalignment, based on sequential Procrustes orthogonal transformations, has been proposed as a method of aligning shared functional information into a common high-dimensional space and thereby improving inter-subject analysis. Though successful, current hyperalignment algorithms have a number of shortcomings, including difficulties interpreting the transformations, a lack of uniqueness of the procedure, and difficulties performing whole-brain analysis. To resolve these issues, we propose the ProMises (Procrustes von Mises-Fisher) model. We reformulate functional alignment as a statistical model and impose a prior distribution on the orthogonal parameters (the von Mises-Fisher distribution). This allows for the embedding of anatomical information into the estimation procedure by penalizing the contribution of spatially distant voxels when creating the shared functional high-dimensional space. Importantly, the transformations, aligned images, and related results are all unique. In addition, the proposed method allows for efficient whole-brain functional alignment. In simulations and application to data from four fMRI studies we find that ProMises improves inter-subject classification in terms of between-subject accuracy and interpretability compared to standard hyperalignment algorithms.

A machine learning based approach towards high-dimensional mediation analysis.

Mediation analysis is used to investigate the role of intermediate variables (mediators) that lie in the path between an exposure and an outcome variable. While significant research has focused on developing methods for assessing the influence of mediators on the exposure-outcome relationship, current approaches do not easily extend to settings where the mediator is high-dimensional. These situations are becoming increasingly common with the rapid increase of new applications measuring massive numbers of variables, including brain imaging, genomics, and metabolomics. In this work, we introduce a novel machine learning based method for identifying high dimensional mediators. The proposed algorithm iterates between using a machine learning model to map the high-dimensional mediators onto a lower-dimensional space, and using the predicted values as input in a standard three-variable mediation model. Hence, the machine learning model is trained to maximize the likelihood of the mediation model. Importantly, the proposed algorithm is agnostic to the machine learning model that is used, providing significant flexibility in the types of situations where it can be used. We illustrate the proposed methodology using data from two functional Magnetic Resonance Imaging (fMRI) studies. First, using data from a task-based fMRI study of thermal pain, we combine the proposed algorithm with a deep learning model to detect distributed, network-level brain patterns mediating the relationship between stimulus intensity (temperature) and reported pain at the single trial level. Second, using resting-state fMRI data from the Human Connectome Project, we combine the proposed algorithm with a connectome-based predictive modeling approach to determine brain functional connectivity measures that mediate the relationship between fluid intelligence and working memory accuracy. In both cases, our multivariate mediation model links exposure variables (thermal pain or fluid intelligence), high dimensional brain measures (single-trial brain activation maps or resting-state brain connectivity) and behavioral outcomes (pain report or working memory accuracy) into a single unified model. Using the proposed approach, we are able to identify brain-based measures that simultaneously encode the exposure variable and correlate with the behavioral outcome.

Differences in functional connectivity distribution after transcranial direct-current stimulation: A connectivity density point of view.

In this manuscript, we consider the problem of relating functional connectivity measurements viewed as statistical distributions to outcomes. We demonstrate the utility of using the distribution of connectivity on a study of resting-state functional magnetic resonance imaging association with an intervention. The method uses the estimated density of connectivity between nodes of interest as a functional covariate. Moreover, we demonstrate the utility of the procedure in an instance where connectivity is naturally considered an outcome by reversing the predictor/response relationship using case/control methodology. The method utilizes the density quantile, the density evaluated at empirical quantiles, instead of the empirical density directly. This improved the performance of the method by highlighting tail behavior, though we emphasize that by being flexible and non-parametric, the technique can detect effects related to the central portion of the density. To demonstrate the method in an application, we consider 47 primary progressive aphasia patients with various levels of language abilities. These patients were randomly assigned to two treatment arms, transcranial direct-current stimulation and language therapy versus sham (language therapy only), in a clinical trial. We use the method to analyze the effect of direct stimulation on functional connectivity. As such, we estimate the density of correlations among the regions of interest and study the difference in the density post-intervention between treatment arms. We discover that it is the tail of the density, rather than the mean or lower order moments of the distribution, that demonstrates a significant impact in the classification. The new approach has several benefits. Among them, it drastically reduces the number of multiple comparisons compared with edge-wise analysis. In addition, it allows for the investigation of the impact of functional connectivity on the outcomes where the connectivity is not geometrically localized.

Sequential fear generalization and network connectivity in trauma exposed humans with and without psychopathology.

While impaired fear generalization is known to underlie a wide range of psychopathology, the extent to which exposure to trauma by itself results in deficient fear generalization and its neural abnormalities is yet to be studied. Similarly, the neural function of intact fear generalization in people who endured trauma and did not develop significant psychopathology is yet to be characterized. Here, we utilize a generalization fMRI task, and a network connectivity approach to clarify putative behavioral and neural markers of trauma and resilience. The generalization task enables longitudinal assessments of threat discrimination learning. Trauma-exposed participants (TE; N = 62), compared to healthy controls (HC; N = 26), show lower activity reduction in salience network (SN) and right executive control network (RECN) across the two sequential generalization stages, and worse discrimination learning in SN measured by linear deviation scores (LDS). Comparison of resilient, trauma-exposed healthy control participants (TEHC; N = 31), trauma exposed individuals presenting with psychopathology (TEPG; N = 31), and HC, reveals a resilience signature of network connectivity differences in the RECN during generalization learning measured by LDS. These findings may indicate a trauma exposure phenotype that has the potential to advance the development of innovative treatments by targeting and engaging specific neural dysfunction among trauma-exposed individuals, across different psychopathologies.

Mode decomposition-based time-varying phase synchronization for fMRI.

Recently, there has been significant interest in measuring time-varying functional connectivity (TVC) between different brain regions using resting-state functional magnetic resonance imaging (rs-fMRI) data. One way to assess the relationship between signals from different brain regions is to measure their phase synchronization (PS) across time. However, this requires the a priori choice of type and cut-off frequencies for the bandpass filter needed to perform the analysis. Here we explore alternative approaches based on the use of various mode decomposition (MD) techniques that provide a more data driven solution to this issue. These techniques allow for the data driven decomposition of signals jointly into narrow-band components at different frequencies, thus fulfilling the requirements needed to measure PS. We explore several variants of MD, including empirical mode decomposition (EMD), bivariate EMD (BEMD), noise-assisted multivariate EMD (na-MEMD), and introduce the use of multivariate variational mode decomposition (MVMD) in the context of estimating time-varying PS. We contrast the approaches using a series of simulations and application to rs-fMRI data. Our results show that MVMD outperforms other evaluated MD approaches, and further suggests that this approach can be used as a tool to reliably investigate time-varying PS in rs-fMRI data.